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^0 \ We present closed analytical expressions for the particle and kinetic energy spatial densities at 

finite temperatures for a system of noninteracting fermions (bosons) trapped in a d-dimensional 
harmonic oscillator potential. For d = 2 and 3, exact expressions for the A''-particle densities are 
used to calculate perturbatively the temperature dependence of the splittings of the energy levels in 
a given shell due to a very weak interparticle interaction in a dilute Fermi gas. In two dimensions, 
• we obtain analytically the surprising result that the Z— degeneracy in a harmonic oscillator shell is 

O ' not lifted in the lowest order even when the exact, rather than the Thomas- Fermi expression for the 

particle density is used. We also demonstrate rigorously (in two dimensions) the reduction of the 
exact zero-temperature fermionic expressions to the Thomas-Fermi form in the large- A*' limit. 
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^ ' I. INTRODUCTION 

a, 

The observation of Bose-Einstein condensation in ultra- low temperature trapped atomic gases (see e.g., Ref. for 
a comprehensive review) , along with the recent experimental effort of Jin and DeMarco |^ in achieving a degenerate 
Fermi gas, has sparked a renewed interest in the study of the thermodynamic properties of trapped quantum gases. 
^ For a dilute gas of trapped bosons or fermions, assuming the interparticle interactions to be very weak, the many- 
body system may be modeled in the zeroth-order approximation as a system of noninteracting particles. This is 
an attractive situation from a theoretical point of view because the thermodynamics of the noninteracting TV-body 
^ ' system can be calculated from the single-particle partition function of the associated trapping potential. For the case 
■ of Bose statistics, exact closed form expressions for the thermodynamic quantities of the ideal trapped gas can be 
\^ obtained at all temperatures in any dimension provided the trapping potential is an harmonic oscillator (HO). An 
, ideal trapped gas of fermions on the other hand is not as easy to analyze. Unlike the Bose gas, the Fermi distribution 
On ' function cannot be expanded as a power series for all temperatures, and as a result, closed form expressions for e.g., 
" the particle and kinetic energy densities, at finite temperature are difficult to obtain. Indeed, it is only recently that 
Brack and van Zyl have constructed exact, analytical expressions for the particle and kinetic energy densities of a 
harmonically confined ideal Fermi gas in any spatial dimension. These expressions are useful for both numerical and 
. analytical computations, but are restricted to zero temperature. Prior to the work in Ref. j^, fermions in harmonic 
' traps have been studied using exact numerical ^, ||, |[ , analytical ^, |l^ and scmiclassical techniques |0 . The 
. exact analytical results for the particle and kinetic energy densities, however, are limited to one dimension (ID) at 
'"^ ' zero temperature or arbitrary dimensions in the high temperature regime p^ . 
Ch 111 this article, our main emphasis will be on presenting exact analytical results for an ideal c?-dimensional trapped 

O Fermi gas at any finite temperature, thereby extending the work in Ref. |^ beyond T = 0. We shall also digress to give 
some novel analogous expressions for the ideal trapped Bose gas which have not yet appeared in the literature. The 
trapping potential will be taken to be an isotropic HO, although our results are easily generalized to the anisotropic 
case. Apart from the intrinsic merit of exact results for a many-body system, our findings should also prove useful in 
^ , the density-functional theory of weakly interacting inhomogeneous Fermi systems. 

' The rest of our paper is organized as follows: In Sec. II, we derive analytical expressions for the finite-temperature 
particle and kinetic energy densities for an ideal trapped Fermi gas in arbitrary dimensions. We also show rigorously 
how the exact quantum expressions at T = asymptotically approach their Thomas-Fermi (TF) results in the large- 
limit; we know of only two other examples where such a result has been rigorously established analytically p^ . 
Then, in Sec. Ill, we derive expressions analogous to the fermionic system for a trapped gas of noninteracting bosons. 
In Sec. IV, we apply our exact result for the spatial density to calculate the mean-field potential, and the resulting finite- 
temperature splitting of the energy levels in a d = 3 HO shell perturbatively. Although a similar zero temperature 
calculation has already been considered numerically by Heiselberg and Mottelson jlj], our results (in contrast to 
Ref. iQ) do not rely on the TF approximation for the particle density. The corresponding calculation for d = 2 yields 
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an interesting result, which is to be discussed in Sec. IV B. Finally, in Sec. V, we summarize our results and suggest 
other possible avenues for future research. 

II. IDEAL FERMI GAS AT FINITE TEMPERATURE 

A. Single-particle density in d dimensions 

The starting point in our investigation is a system of nonintcracting fcrmions at zero temperature described by the 
time-independent Schrodinger equation 

H(t>,{r) = [f + y(r)]4(r) = e,0,(r) , (1) 

where ^(r) is a one-body potential to be specified later (all e^'s are taken to be positive). The single-particle density 
matrix can be obtained by an inverse Laplace transform of the zero-temperature Bloch density matrix, Co(r,r'): 



|co(r,r';/3) 



(2) 



p(r,r') = 2 ^ 0*(r')0.(r)e(i?/) = 

e,<Ef 

where 

Co(r,r';/3) = ^ exp(-/3£,) , (3) 

all i 

and Ef is the Fermi energy; the factor of two accounts for spin. We have put in the unit step-function Q{Ef) in 
Eq. (H) so that the Laplace transform with respect to Ef may formally be taken to be two-sided Note that in 
quantum statistical mechanics, P is usually identified with the inverse temperature, P = l/fc^T. However, in our 
present context, P is to be interpreted as mathematical variable which in general is taken to be complex, and not the 
inverse temperature l/fcsT. 

At finite-temperature, the single-particle density matrix is obtained from the Bloch density matrix by using the 
relation ml 



p(r,r';T) = /:, 

where 



|c^T(r,r';/3) 



(4) 



CT(r,r';/3) = Co(r,r';/3)-^^ , (5) 

sm(7r/5i ) 

is the finite-temperature Bloch density matrix, and fi is the chemical potential. In Eq. (4) , the Laplace transform with 
respect to fj, is two-sided, so that fi is allowed to go negative. Specializing now to the case of an isotropic harmonic 
oscillator in d dimensions, viz., 

V{r) - ^niLuy, r - ^ xj + + ■ ■ ■ + x} , (6) 
we have for the zero-temperature Bloch density matrix 

C^^)(r,r';/3) = c('^)(g,s;/3) = (^^y^' exp {- [q^ tanh(/3/2) + (5V4) coth(/3/2)] } . (7) 



In the above expression (and what follows), all lengths and energies have been scaled by Iq = a/KT™^ ^^'^ ^wq, 
respectively, and we have introduced the center-of-mass and relative coordinates: 

q=i(r + r'), s = r - r' . (8) 

In order to obtain the finite-temperature single-particle density, we need to set s = in Eq. (^ and along with 
Eq. (||), perform the inverse Laplace transform given by Eq. (|). This is a difficult task if one attempts to use Eq. (m) 
in its present form. However, matters can be simplified considerably if one uses the following identity |8j, 

00 

exp{-:Etanh(/3/2)} = ^(-l)"L„(2a;)e-"{e-"'^ + e-("+i)'3} ^ (9) 

ra=0 
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in Eq. which for s = now reads 

Substituting Eq. ( [lo[ ) into Eq. (||) and performing the inverse Laplace transform, Eq. (^), leads to the finite- 
temperature single-particle density in any dimension. In order to illustrate the method, we will now proceed to 
give an explicit calculation for the case of two dimensions (2D), followed by a statement of the general result in 
arbitrary dimensions. 



1. Two dimensions 



The 2D system turns out to be the simplest case of all, and for this reason, we give here a detailed derivation 
of the finite-temperature single-particle density. We begin by noting the following important exact inverse Laplace 
transforms (all two-sided): 



'-"q 



-nl3 



sinh(/3) 



fc=0 



ttT 



sin(7r/3r) 



[exp(-f) + l] 



Putting d = 2 in Eq. ( [lO| ) and using Eqs. (||), (||), the finite-temperature density is given by 



ri=0 



-n/3 _^ g-(n+l)/3\ 



sinh(/3) J sm{TTf3T) 



(11) 
(12) 

(13) 



Applying the convolution theorem for Laplace transforms ||l^ and making use of Eqs. (|T^), (p^), we immediately 
obtain 



p(2)(q;r) 



^(-l)"i„(2(72)e-'?' X ^ ] / drSir - {2k + 1) 



1 



[exp(l^) + 1] 



+ / dT6{T - {2k + 2)-n) 



[cxp(l^) + 1] 



= ( I) V(-l)"i„(2g2)e-<?^ X V J J ^ + J 

V-/,1^o .t^\[exp(2MlLtWi) + i] [exp(H^±!^^ 



5](-l)"L„(2g2)e-^^ x ^ 



[exp(^ 



where the function Fi^-* (/i) is defined as 



^0 [exp(S-^^) + 1 



(14) 



(15) 



and £„ = {n + 1) is the noninteracting energy spectrum (in scaled units) of the 2D harmonic oscillator potential. In 
the limit T — > 0, the Fermi function in Eq. (|l5| ) goes over to the Heaviside step function, and when filling M + 1 
oscillator shells, Eq. (|lj) reduces to 



M 



P^^Hq)=(-)Y.(M-n+l){~irLU2q^)e-^ 



(16) 
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FIG. 1: The 2D exact (solid curves) and TF (dashed curves) number densities at various temperatures for A'' = 420 particles. 
The upper curves in the main figure correspond to T = 0, while the lower-lying curves are for T = 5. The figure inset shows 
the 2D exact and TF densities at T = 1. All quantities have been scaled as discussed in the text. 



Equation ( p^ ) is of course identical to the result obtained in Ref. |^ in which only zero-temperature quantities were 
investigated. In Fig. |l|, we display the exact 2D radial density (solid curves), p^^^(r; T), for N — 420, where the lowest 
20 shells are filled. For comparison, we have also displayed the corresponding 2D TF densities (dashed curves) at zero 
and finite temperatures We have used Eqs. ( |l^ and (^6|) for the exact nonzero and zero-temperature density 
respectively. In all plots, T denotes the dimensionless quantity (^^j-). We note from Fig. |^ that the oscillations due 
to shell effects are already washed out by T = 1, resulting in TF and exact densities being almost indistinguishable. 
In the high temperature regime where ^ < Sq, we can make use of the expansion 



exp( 



T 



3 + -^ ^3t^/T -3(n+l)/T ~2jk/T 



(17) 



J = l 



in Eq. (^^. The k sums can be evaluated exactly, and we obtain 

•,(2) 



p^-'(q;T) = 2^(-l)^+ie-'"^/^Cf (q;j7r) 



(18) 



Making the formal identification oi (3 — 1/T, Eq. (081) can be re- written as 



p(2)(q;T) = 2^(-l)^"+ie^"'^^cf (q;j/3) 
i=i 



(19) 
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A result identical to Eq. has recently been obtained by Wang pO[ , but as we have discussed above, this form for 
the density is only valid for < Eq. In contrast, our Eq. (|l4| ) is valid at all temperatures and appears in a form that 
clearly displays the role of the quantum statistics of the system. 

Finally, an interesting relationship between Fn (/i) and the Fermi distribution function can be obtained through 
the normalization condition 

N = I d2qp(2)(q;T) 

OG 

^ (20) 

,£^Jexp(^) + 1] 



2. Arbitrary dimensions 

The generalization of the above analysis to arbitrary dimensions is straightforward. First, since we are only inter- 
ested in the diagonal Bloch density matrix, we can immediately set s = in Eq. (j^). Furthermore, we observe that 
in any dimension 

sinh(i-^/2)(/3) = 2W2-i)e'5(i-^/2) (l+f^ sL'^e'^^'^^ . (21) 

\ m=l / 

Multiplying the numerator and denominator of Eq. by Eq. (^l|), and using the identity, Eq. (^), we obtain 



n=Q 

Proceeding exactly as in the 2D case, we finally obtain in any dimension 

d/2 OO 

TT 



p('^)(q;T) = 2(-) ^ FW(A.)(-l)"L„(2g2)e-^\ (23) 

where 



n=0 



^ Vexp[(er +k- fi)/T] + 1 exp[(e;"^ +k + 2m^ fi)/T] + 1 ^ 



and eif^ = n + d/2. The expansion coefficients for c? = 1, 2, 3 are given by H 

gl'J = -(2m - 3)!!/(2m)!! with g^"^ = -1/2 , 
g^n^ = , 

gi'J = (2m-l)!!/(2™)!! . (25) 

One of the most appealing aspects of Eq. ( p3| ) is that the function Fii'\^) contains all of the dimensional and 
statistical dependence. Moreover, for temperature and particle number ranges where significant deviations from the 
TF approximation are expected (see Fig. |^), the above sums can be truncated quite quickly. For example, at low 
temperature and N — 0(10'^), (rimax,kmaxT'mmax) = (50,50,10) is more than adequate to obtain single-particle 
density in any dimension. 

Following the 2D example, it is easy to show that in any dimension, the high-temperature density (i.e., /i < Eq''^) is 
given by 

OO 

pW(q;T) = 25](-l)^"+ie^-'5^c('')(q;i/?) . (26) 
i=i 
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In fact, when fi < £'^\ the following identity for fermions exists in all dimensions: 



[2 sinh(j/3)]'i/2 

In the d-dimensional case, the chemical potential is determined through the normalization condition 

TV^"*) = J d'^'^\p[q,T) 

OO 

= 2^(-irf^W(/.)xFhn,d/2;l;2] 



(27) 



2E 



(28) 



where F[a,b]c; z] is a hypergeometric function p^ , and ry^''' denotes the degeneracy of the quantum levels in d 
dimensions. The functions -F[a, 6; c; z\ can only be evaluated analytically for even d; their values are given by 
(-l)"F[-n,d/2;l;2] = 1, (2n+ 1), (2n2 + 2n+ 1), for (i = 2,4,6, respectively. 

B. Single-particle kinetic energy density in d dimensions 

The finite-temperature kinetic energy density can be obtained in a manner entirely analogous to the derivation for 
the single-particle density. We begin by first investigating two expressions at zero temperature (in scaled units) 



T(q) = -2i '^»WVV.(r) , 



(29) 



ri(q)=2l I^<^^«l' ' 



(30) 



ei<Ef 



where, again, the factor of two accounts for spin. Equations (|2^), (^0|) both integrate to the same exact kinetic energy 
and in the presence of time reversal symmetry are simply related by 



T(q) =ri(q)- ivV(q) 



(31) 



At low temperatures, shell effects give rise to oscillations in r(q) and ri(q) that are exactly opposite in phase [gO| 
As such, a convenient quantity to consider is their mean: 



^(q) = l^[r{^)+rM)] 
= -^V2p(q,s)|,=o 



(32) 



Starting from Eq. ^ and noting that differentiation commutes with the Laplace inverse operator, we readily obtain 
after suitable manipulations of hyperbolic functions the following Laplace inverse: 



4'^(q;/3) 



coth(/3/2) 



(33) 



where c'^\c\;, (3) is defined by Eqs. (||, |^), with s = 0. This immediately gives 



e(^)(q;r) 



2 M 



4sinh^(/3/2) 



.sinh(l-'^/2)(^)g-9=tanh(;3/2) 



sin(7r/3r) 



(34) 
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In order to proceed, we su bstit ute Eq. into ( |34| ) and make use of the convolution theorem for (two-sided) Laplace 
transforms just as in Sec. [I A. We require the following exact Laplace inverse 



-n(i 



4sinh^(/3/2) 



= ^(fc + l)(5(r; - (2fc + 1) - n) 6(7?). 



fe=0 



A direct evaluation of Eq. ( J34| ) leads to the simple expression for ^(q; T), which is valid in any dimension: 

1 d 



^G('^)(M)(-irL„(2g2), 



n=0 



where 



fc=i 



exp[(el'^^ + fc - 1 - ^)/r] + 1 ^1 exp[(e5°^ + + 2m - 1 - /i)/T] + 1 



E 



(35) 



(36) 



(37) 



It is interesting to note that aside from numerical factors, the finite-temperature kinetic energy density and single- 
particle density have identical functional forms in all dimensions. 

An alternative, and useful formula for the kinetic energy density, ^^'*)(q;T), may be obtained from Eq. ( ^ by 
noting that 



g-,«ta„h(^/2)^/^g/ = icoth(/3/2) e-9''-^h(;3/2)_ 

It then follows that the kinetic energy density is given by 

C(^)(q;r) = -y p<^'\ci';T)q'dq' . 



(38) 



(39) 



Reduction to the TF form for N ^ 00 



At low temperatures (i.e., T = 0) and small particle numbers, the spatial densities, given by Eqs. (p3|,36), exhibit 
pronounced shell oscillations, especially in low dimensions However, as the number of particles is increased, these 
shell effects get washed out and the spatial densities approach their usual TF form. In this section we show rigorously 
how the exact T — density profiles reduce to their TF forms as iV — > 00. For simplicity, we focus only on the 2D 
case. 

Consider the exact zero-temperature density given by Eq. (|l6|). Putting 2q^ — t, and using the one-sided Laplace 
transform 



(40) 



it is easy to find that 



M 



Fis) 



Tl = 



(i + s)»+i 



(41) 



The above n sum may be performed exactly pT[ , yielding 



Fis) 



2 1 



1 - 2s 
1 + 2s 



M+l 



{l-2s} -I- (4sAf -l-6s- 1) 



(42) 



2s-l 



2s+l 



< 1, 



Note that the real part of s must be positive when taking the inverse Laplace transform of F{s). Hence 

so that the first term in Eq. ( ^2| ) vanishes as M ^ 00. This is indeed the case when the number of fermions N ^ 00. 
In this limit, we get 



Fis) 



1 



M 



TTS 2tis Atts"^ 



(43) 
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Taking the inverse Laplace transform with respect to t, we immediately obtain 

The last expression is the large-N asymptotic TF result for the density in d = 2 dimensions. Its normalization to the 
correct particle number TV = + 3M + 2 for M + 1 filled shells leads to the Fermi energy Ep = M + 3/2 + 0(1/M) 
in the large-M limit. 

An analogous calculation may be performed for the kinetic energy density at T = 0. Taking d = 2 in Eq. (^), and 
going over to the T — limit, we have 

1 ^ 

^(2)(q;T = 0) = - V(-l)"(M-n+l)2L„(i)e-*/2 . (45) 

TT ^ — ' 

n=0 

As before, we take the one-sided Laplace transform of the above equation with respect to s, and perform the n sum 
exactly. We then take the large M limit, Laplace invert back to the t = 2q^ variable and obtain (M + 3/2 — (7^/2)^/27r, 
which is just the TF result for the kinetic energy density. 



III. IDEAL BOSE GAS AT FINITE TEMPERATURE 



A powerful aspect of formalism can be illustrated when one considers a trapped ideal Bose gas at finite temperature. 
In this case, the finite-temperature Bloch density matrix is simply given by (compare with Eq. (5)) 



CT(q,s;/3) --Co(q,s;/3) 



TTf3T 



tan(7r/3r) 



Noting the following inverse Laplace transform 



nT 



tan(7r/3T) 



[exp(-^)-l] 



(46) 



(47) 



we observe that Bose statistics serve only to change the behaviour of the functions F^\^) and Gif \^j,). Explicitly, 
for the Bose gas we find 



/ -,\d./2 oc 

(q;r)= - ^FW(M)(-l)"L„(2g>-« , 



n=0 



where 



Vexp[(£i'*^ +k- n}/T] - 1 exp[(e}r' +k + 2m- n)/T] - 1 



E 



and similarly. 



(48) 



(49) 



1 d 



^^'Hci;T) = _-^G('^)(^)(-l)"L„(2g2)£ 



(50) 



where 



k=l 



exp[(ei''^ + fc - 1 - fi)/T] - 1 ,,^1 exp[(e^"^ +k + 2m~l~ fi)/T] - 1 



ym 



(d) 



(51) 



Just an in Sec. all of the dimensional and statistical dependence is buried in the functions F^''-' and Gif\ Not 
surprisingly, the Bose versions of F^\ii) and Glf\ii) differ from their Fermi counterparts only by the sign in front 
of unity in the denominator. 
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1. Alternative method 

As it happens, there is an alternative method for evaluating the finite-temperature spatial density in the case of 
Bose statistics which avoids the summations introduced by the inverse Laplace method. One starts by noting that at 
all temperatures, the Bose distribution function can be expanded as a power series 

oo 

This immediately allows us to write the finite-temperature density matrix in the form (here, the bosons are taken to 
be spinless) 

oo oo 

p('^)(r,r';T) = ^ e^^T ^ (r')</.„(r) exp(-M'VT) ■ (53) 

j=l n=0 

Identifying (3 — 1/T in Eq. (3) then allows us to write Eq. (53) as 

oo 

(q, s; T) = ^ e^P^^Ci'^ (q, s; jP) . (54) 

The finite-temperature particle density is obtained by setting s = in Eq. (54) and the kinetic energy density is once 
again evaluated from Eq. (32). 

It is straightforward to show that our expressions for the Bose gas obtained by the inverse Laplace method lead to 
results identical to the so-called alternative method. For example, in the particularly simple case of two dimensions, 
we have 



/ 1 \ °° oo f 

^"^^ n=0 fc=0 [ [e 



exp( '^+"+^-^ ) - 1 exp( "^+"+"-^ ) - 1 



r2k+n+2-fj,) 



(55) 



Expanding the Bose functions as in Eq. (52) and performing the k sums exactly, we obtain {(3 = 1/T) 



P 



(2) 



oo 

= Y.e=P^CM.3P), (56) 

by virtue of Eq. (9). Indeed, in all dimensions and all temperatures, the following identity can be established for 
bosons 



[2 sinh(j/3)]'i/2 

where F^' {p) is given by Eq. (49). The analogous result for fermions is given by Eq. (27) 



IV. WEAKLY INTERACTING FERMI GAS 

Up to this point, we have been considering noninteracting gases in a trap. In this section we shall investigate 
(following Ref. iQ) a very dilute weakly interacting Fermi gas a finite temperature. We will consider both three- 
dimensional (3D) and 2D geometries and examine the effect of the interaction on the degeneracy of levels in a given 
HO shell. 
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FIG. 2: Finite-temperature energy shift, Eq. (y), for a weakly interacting 3D trapped Fermi gas as a function of the scaled 
temperature T. The sohd curves correspond to np — 10 and the dashed curves to np ~ 50. Note the weak temperature 
dependence on the energy shift for np = 50. 



A. Three dimensions 



We first consider iV fermions in a 3D isotropic harmonic oscillator trap interacting via a contact pseudo-potential. 
The Hamiltonian is then given by 

^ / 2 1 \ 

^ - E Y + + '™ ^ ' • ^^^^ 

i=l ^ ^ i<j 

In the above, a is the dimensionless s-wave scattering length. The mean-field potential at temperature T is 

U{r;T) = 2Trap{r;T), (59) 

where, keeping in mind that d = 3, we omit the superscripts on all quantities. If the system is very dilute in the sense 
that hiv » \U\, we may use the unperturbed expressions, Eq. (p3[), for p{r;T) to obtain U{r;T). To fix the notation 
for the degenerate states in a shell, we consider the noninteracting system at T = 0, for which the Fermi energy is 

EF^inF + i/'i) , (60) 

and Tii? is the principal HO quantum number at the Fermi surface. The quantum number np- is related to the number 
of fermions, iV, in the well by 

N ='^{nF + l){nF + 2){nF + 'i) . (61) 
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FIG. 3: Finite-temperature energy shift, Eq. ( |64[ ), for the 3D weakly interacting Fermi gas as a function of np- The open circles 
correspond to / = np whereas the filled circles are for Z = or Z = 1. For large np, F(nF,l;T) approaches a constant value, as 
is expected in the TF limit. 



Note that the shell specified by np has degenerate states with angular momenta I — np,np — 2,np — 4, 1 or 0. 
Consider now a particular radial HO state TZni, with (2n + I) = np in this shell. Due to the mean-field potential 
U{r;T), the single-particle energy of this state is shifted (in lowest order perturbation theory) by the amount 



Ep = A. 



^^)|7^„^|V 



dr 



(62) 



Thus the degeneracy of the various angular momentum states is lifted by the temperature-dependent interaction 
(p9|). Using Eqs. (|5^) and (^3|), we note that the relevant radial integral to obtain the energy shift is of the form (at 
any temperature T) 



I 



^21+2^- 



dr 



(63) 



Although it is possible to obtain an analytical expression for the above integral, in practice it is more convenient to 
calculate it numerically. In Fig. ||, we display the results for the energy shift as a function of the scaled temperature 
T for I = and I = np, with np = 10 and 50. Since this energy shift obviously depends on the scattering length a 
and np, we choose to plot the dimensionless quantity 



F{np,l;T) 



Ae„^;(r) 
27raiVi/2 



(64) 



Our definition of F{np, I; T) in Eq. ( |64| ) is a little different from that in [Q, where the calculation was performed at 
T = in the TF approximation. Note that the temperature dependence for nj? = 10 is much more pronounced than 



12 



for np = 50. This would be less so if one scaled the temperature in units of the Fermi energy rather than the oscillator 
spacing. We also show in Fig. || the level splittings, F{np,l;T), as a function of np at various temperatures. For 

large np, we observe that F{np,l;T) approaches a constant, demonstrating that Ae„j^z(T) behaves as iV^/^ ~ n^J^ , 
as is to be expected in the TF limit. 



B. Two dimensions 



We now calculate the corresponding energy shift in the quasi-2D harmonic oscillator trap. A highly anisotropic 3D 
oscillator will behave as such when the oscillator frequency in one direction (say ujz) is much, much greater than in the 



other two {ujx 



UJv 



Furthermore, the temperature T under consideration is low enough so that ksT <C hojz 



If the interatomic interaction is so weak that it cannot cause excitations in the z direction, we may restrict the Hilbert 
space in the z direction by setting the oscillator quantum number = 0. The effective interaction of the quasi-2D 
system is then obtained by taking the 3D delta- function potential in Eq. (p8|), and taking its expectation value with 
respect to the HO wave function in the z direction. 



Mz) = (^)'^^ exp(-c^,zV2) . 
The resulting quasi-2D Hamiltonian is then given by (again, scaled units are used) 

^ / 2 1 \ 



(65) 



(66) 



where the momenta and coordinates are planar vectors, g — 2v27r(a/£^) is the effective 2D dimensionless coupling 
constant, and we have denoted the oscillator length in the z direction by Iz- We restrict this calculation to T = 0, 
although the extension to finite temperature is straightforward. The mean-field one-body potential is then given by 



where the 2D p^'^\r) is given by Eq. (|l^). Identifying M with np, we note that for the 2D case, 

N = {np + l){np + 2) , Ep = {np + 1) . 
The radial 2D HO wavefunction is given by 



2nl 



in+\l\)l 



1/2 



with the normalization drr\Rni\'^ — 1. The energy shift analogous to Eq. (62) is 



drrp{r)\R„i 



and 



F'^'HnpJ) 



Ae 



(2) 

np ,1 



(67) 



(68) 



(69) 



(70) 



(71) 



The 2D energy shift at T = is shown in Fig. ^. It readily seen from this figure that F'^'^\np,l) depends on np-, but 
not explicitly on I. Consequently, unlike the 3D case, there is no splitting (in the lowest order perturbation theory) 
between states with different I values in a given shell at T = 0. This surprising result may be also checked analytically. 
Specifically, a little algebra gives 



Ae, 



(2) 



2 n\ 

2a/- g- 



TT ''in+\i\y. 



J2 (-1)'"("F - TO + 1) / da:a;l'le-2^L™(2x) [ij^l (x) 



(72) 
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FIG. 4: The energy splitting F'-'^\nF,l) defined by Eq. (71) for a 2D Fermi gas at T = 0. Note that F^'^\nF,l) is independent 
of I for every HO sheU. 



The integral above can be performed exactly by noting that pi 



{n+\l\)\ 



E 

/c=0 



J2k 



(I^l + fc)!fc! 



(73) 



The final result is given by (with the restrictions {n — k) < m < {n + k + \l\)) 



Ae 



(2) 



— riF 



+ nvHl^ (l^l + 2fc)!(n + fc + |;|)! 

" ^ t:^ 2l'l+2'= kl{n - k)l{\l\ + k)l{m - n + ky.{n + k + \l\ - m)\ 



(74) 



An explicit case by case evaluation of the above expression verifies that the shift depends only on np, and is 
independent of |/|; we have not been able to demonstrate this fact analytically for a general np. As stated earlier, 

the independence of Ae^^^ ^ on \l\ is to be expected for large np in which the TF density depends only on r^. What 
is surprising is that the shift is independent of / for all I within a given shell even though the exact single-particle 
density, Eq. (16), has been used. 



V. SUMMARY AND DISCUSSION 



We have presented for the first time exact analytical expressions for the number and kinetic energy densities of an 
ideal Fermi gas in a harmonic trap at finite temperatures. These new expressions, given by Eq. (^3|) and (|36|), are 
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valid at any temperature, and in any dimension. We have also derived the corresponding results for bosons, which 
are given by Eqs. ( ^8|) and (|50|). UnHke previous expressions in the literature (see e.g., Ref. 10), our functional forms 
for the particle and kinetic energy densities clearly illustrate the role of the quantum statistics and have a universal 
form independent of the dimension of the system. We have also demonstrated analytically in Sect. II C (for d — 2), 
how the exact fermionic expressions for the particle and kinetic energy densities at T = reduce to their TF forms 
as iV — > oo. It is worthwhile noting that establishing this reduction analytically using a wavefunction representation 
for the particle and kinetic energy densities is a non-trivial task (see e.g., Ref. [^). 

To complement our ideal gas results, we have also considered a weakly interacting trapped Fermi gas. The exact 
fermionic expression for the spatial density was used to obtain the mean-field potential, and the finite temperature 
first-order shifts in the energy levels for a given HO shell. Similar calculations were recently reported by Heiselberg 
and Mottelson in 3D, but only at T = and within the TF approximation. As was pointed out in Ref. for 
attractive interparticle interactions (i.e., negative a), there are also pairing correlations between time- reversed states 
(Z,m) and (/, — m) in a given shell. In the situation when the splittings given by Eq. ( |62| ) are small compared with 
the pairing gap, the pairing may take place between all states of a shell leading to supergaps. We plan on performing 
a detailed study of the supergap phenomenon at finite temperatures in a future publication. 

As our final application, we also studied the energy level splittings in a d = 2 isotropic HO geometry at T = 0. We 
found the surprising result that the degeneracy for a given HO shell is not lifted in the lowest order perturbation 
theory, even when the exact expression for the single-particle density is used. In the 2D TF approximation, this result 
can be understood by noting that the TF particle density, p"r^\{r) oc (/i — r^/2), does not break the SU{2) symmetry of 
the harmonic trap. However, the radial dependence for the exact density (especially for a small number of fermions) 
is more complicated, and therefore our result is unexpected. Given this finding, it should prove very interesting to 
study the nature of pairing transitions in a 2D trapped Fermi system at finite temperature. 
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